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INTERCOMPUTER TRANSFER IN FULL PRECISION OF 
ARBITRARY DATA ON MAGNETIC TAPE EMPLOYING 
NASTRAN USER TAPE FORMAT 

By James L. Rogers, Jr. 

Langley Research Center 


SUMMARY 

A description is presented of two new utility programs which implement the trans- 
fer, in full precision, of arbitrary data (matrices or tables) between any of the three 
NASTRAN operative computers without the handling of large card decks. These com- 
puters include the CDC 6000 series, the IBM 360-370 series, and the UNIVAC 1100 series. 
The data may be generated by NASTRAN or by another computer program if the NASTRAN 
user tape format is employed. 


INTRODUCTION 

The NASTRAN computer program (refs. 1 and 2) is capable of execution on three 
different types of computers, namely, the CDC 6000 series, the IBM 360-370 series, and 
the UNIVAC 1100 series. A typical activity requiring transfer of data between dissimilar 
computers is the analysis of a large structure such as the Space Shuttle by substructuring. 
Models of portions of the vehicle which have been analyzed by subcontractors on their 
computers must be integrated into a model of the complete structure by the prime con- 
tractor on his computer. Presently the transfer of NASTRAN matrices or tables between 
two different types of computers is accomplished by punched cards or a magnetic tape 
containing card images. These methods of data transfer do not satisfy the requirements 
for intercomputer data transfer associated with a substructuring activity because (1) accu- 
racy will be lost due to the precision limitations (10 significant digits) of the NASTRAN 
Direct Matrix Input (DMI) punched card and (2) large order matrices make card handling 
too cumbersome. 

To provide a more satisfactory transfer of data, two new programs, RDUSER and 
WRTUSER, were created. These two programs, used in conjunction with NASTRAN mod- 
ules OUTPUT2 (ref. 1, p. 5.3-20h) and INPUTT 2 (ref. 1, p. 5.3-16h) available in level 15 
and later versions of NASTRAN, allow data to be transferred between computers without 
loss of accuracy and without handling large decks of punched cards. The- purpose of this 
paper is to describe both the utility programs RDUSER and WRTUSER and their applica- 



tions by typical data transfer. Although data may come from any computer program 
using the NASTRAN user tape format, examples in this paper are confined to NASTRAN 
data since RDUSER and WRTUSER were written with the NASTRAN user in mind. 

OVERVIEW OF PROGRAMS 

Beginning with level 15 NASTRAN provided the capability of using FORTRAN 
WRITE statements to write intermediate data blocks (matrices or tables) on a magnetic 
tape. This was made possible by the NASTRAN module OUTPUT2 which has the follow - 
’ ing calling sequence: 

OUT PUT 2 DB1, DB2, DB3, DB4, DB5//V, N, Pl/V, N, P2/V, N, P3 $ 

where the DBi are the data blocks to be written on tape, PI is a parameter for positioning 
the tape, P2 is the FORTRAN unit number assigned to the tape, and P3 is the FORTRAN 
User Tape Label (default = XXXXXXXX). 

The tape created by OUTPUT2 is a binary tape. Its format is shown in tables I, II, 
and III. Tapes created by a program other than NASTRAN are acceptable as long as the 
data are output in this format. In order to write the header information on the tape, the 
PI parameter must be -1 (rewind before writing) the first time OUTPUT2 is called in 
NASTRAN; otherwise, PI is zero. This binary tape must be converted to a BCD (binary 
coded decimal) tape before it can be used on a computer of a different type. The conver- 
sion is performed by the utility program RDUSER which accepts tables and single- or 
double-precision, real or complex matrices. No precision is lost in generating the BCD 
tape, and the problem of handling large numbers of punched cards is alleviated. 

The tape containing the BCD data is transferred to another installation. Before 
these data can be used as input for NASTRAN at this installation, two tasks must be per- 
formed. The first task is to convert the source of the BCD tape, written by RDUSER, to 
another source form readable by the computer on which the data will be used. The sec- 
ond task is to convert the BCD tape into an acceptable binary form for the NASTRAN mod- 
ule INPUTT2. The program WRTUSER accomplishes both of these tasks. The calling 
sequence for the INPUTT2 module has the form 

INPUTT 2 DBI, DB2, DB3, DB4, DB5/V, N, Pl/V, N, P2/V, N, P3 $ 

where the DBi are the data blocks to be recovered from the binary tape, PI is a parameter 
for positioning the tape (PI must be -1 for the first call to INPUTT2 and zero on all suc- 
ceeding calls), P2 is the FORTRAN unit number assigned to the binary tape, and P3 is the 
FORTRAN User Tape Label (default = XXXXXXXX). A flow chart of the complete tape 
interface method is as follows: 
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DESCRIPTION OF RDUSER AND ITS SUBPROGRAMS 


Program RDUSER 

RDUSER is a FORTRAN main program that has as its primary function the calling 
of various subprograms for data manipulation. The flow chart and the program listing 
that follow show how RDUSER controls program operation. 
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o n o n n n o n o o non non n n n non non non non n n n n 


THIS PROGRAM CONVERTS A BINARY UNFORMATTED TAPE OUTPUT FROM 
NASTRAN MODULE OUTPUT2 INTO A BCD TAPE 

DIMENSION I RD I N (7 ) 

DATA I BLANK, I STOP. I T/4H ,4HST0P,4H IT / 

WRITE (6,160) 

COMENT READS USER INFORMATION ABOUT DATA BLOCKS ON TAPE 
CALL rdcom 
read in tape id 


10 READ (3.90) NAM 1 , NAM2 

IF (NAM1 .EQ» IBLANK.AnD«NAM2,EQ, IBLANK) GO TO 80 
I SW = 0 

READ tape until MATCHING ID IS FOUND 
20 READ <9) I CNT 

IF ( I SW.EQ. 0 » AND. I CNT »NE .3 ) GO TO 70 
ISW=1 

READ (9) ( IRDIN< I ), 1=1 , I CNT) 


WRITE INFORMATION ON 
IF ( I CNT ,EQ, 3 ) WRITE 
WRITE INFORMATION ON 
IF (ICNT*E0,7) WRITE 
WRITE INFORMATION ON 
IF (ICNT.EQ.2) WRITE 


BCD TAPE 

(8,100) ( I RD I N ( 1 > « I x 1 * I CNT ) 

BCD TAPE 

(8,110) ( IRD!N( 1 ) , 1=1 , ICNT) 

BCD TAPE 

(8,120) ( IRDIN( I ) , I =1 , ICNT) 


CHECK FOR MATCH 

IF (NAM1 »EQ, I RD IN ( 1 ) • AND,NAM2»EQ, I RD IN ( 2 ) ) GO TO 50 
GO TO 20 

50 READ (9) IDUM 
READ (9) IDUM 

ID HAS BEEN FOUND ON TAPE 
I CODE® 1 WE ARE LOOKING FOR A MATRIX 
IC00E=2 WE ARE LOOKING FOR A TABLE 
IOPT=0 DO NOT PRINT DATA BLOCK ELEMENTS 
I OPT* 1 PRINT DATA BLOCK ELEMENTS 

60 READ ( 5,90) NAM1 « NAM2 . I CODE , I OPT 

IF (NAM1 ,EQ, IBLANK, AND.NAM2.EQ, IBLANK ) GO TO 10 
IF ( I CODE .EO • 1 ) CALL MATRD ( NAM 1 , NAM2 , I ERR , I OPT ) 

IF ( 1 CODE.EQ .2 ) CALL TABRD (NAM 1 , NAM2 . I ERR , I OPT ) 

IF (IERR.EQ.O) GO TO 60 

MATRIX OR TABLE NAME NOT FOUND ON TAPE 

WRITE (6,130) NAM 1 , NAM2 
GO TO 80 


RDU 1 1 
RDU 12 
RDU 13 
RDU 14 
RDU 20 
RDU 30 
RDU 40 
RDU 50 
RDU 60 
RDU 70 
RDU 80 
RDU 90 
RDU 100 
RDU 110 
RDU 120 
RDU 130 
RDU 131 
RDU 140 
RDU 150 
RDU 160 
RDU 170 
RDU 180 
RDU 181 
RDU 190 
RDU 210 
RDU 220 
RDU 230 
RDU 240 
RDU 250 
RDU 260 
RDU 270 
RDU 280 
RDU 290 
RDU 300 
RDU 310 
RDU 320 
RDU 330 
RDU 340 
RDU 350 
RDU 360 
ROU 370 
RDU 380 
RDU 390 
RDU 400 
RDU 410 
RDU 420 
RDU 430 
RDU 440 
RDU 450 
RDU 460 
RDU 470 
RDU 480 
RDU 490 
RDU 500 
RDU 510 
RDU 520 
RDU 530 
RDU 540 
RDU 550 
RDU 560 
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non 


tape id was not found on tape 

70 WRITE (.6*140.) NAM 1 *NAM2 
STOP * 

C 

C WRITE EOF ON TAPE 

C 

80 I CNT=0 

WRITE (.6*1 50 ) ICNT 
WRITE (8*170) I STOP * IT 
WRITE (6*160) 

STOP 

C 

90 FORMAT <2A4*2I2) 

100 FORMAT (SOX* 3110) 

110 FORMAT ( 1 OX * 7 ( A4 * 6X ) ) 

120 FORMAT ( 6 OX « 2 ( A4 * 6X ) ) 

130 FORMAT (1 OH THE NAME *2A4 *23H WAS 
140 FORMAT ( 7H LABEL *2A4 * 23H WAS NOT 
150 FORMAT (70X*I10) 

160 FORMAT (I HI) 

170 FORMAT (72X*2A4) 

END 


ROU 570 
RDU 580 
RDU 590 
RDU 600 
RDU 601 
RDU 607 
RDU 608 
RDU 609 
RDU 610 
RDU 620 
RDU 621 
RDU 630 
RDU 640 
RDU 650 
RDU 660 
RDU 670 
RDU 680 
RDU 690 

NOT FOUND ON TAPE /) RDU 7q0 

FOUND ON TAPE /) RDU 710 

RDU 720 
RDU 730 
RDU 731 
RDU 740- 
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Subprogram RDCOM 


RDCOM is a FORTRAN subprogram. The primary job of RDCOM is to read com- 
ments by the user about the data on the tape. These comments are punched on cards for 
input into RDUSER. The cards are read with a free field format and thus may contain 
any information describing the data blocks. The comments are written On unit 8. The 
RDCOM subprogram listing follows. 



SUBROUTINE RDCOM 

THIS SUBROUTINE READS COMMENTS FROM 

DIMENSION I COM (20) 

DATA 1QU1T/4H END/ 

READ (5.20) (ICOM< I ) ,1=1 *20) 

PRINT 20* ( I COM < I ) ' 1 = 1 *20 > 

WRITE ( 8 . 20 ) < ICOM( I ) . I* 1 .20 ) 

IF ( ICOM (20 ) .NE» I QU1 T ) GO TO 10 
RETURN 

FORMAT (20A4) 

END 


COM 1 0 
COM 1 1 

CARDS AND WRITES THEM ON TAPE COM 12 

COM 1 3 
COM 20 
COM 30 
COM 40 
COM 50 
COM 60 
COM 70 
COM 80 
COM 90 

COM 110- 
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Subprogram MATRD 

MATRD is a FORTRAN subprogram. Its parameters are the matrix name (NAM1, 
NAM2), an error flag (IERR), and a print option (IOPT) for the user to print or not to 
print the elements of a column. The primary job of MATRD is to read the matrix data 
from unit 9, unpack the columns, and write them on unit 8 in a suitable format. If the 
matrix name cannot be found on unit 9, IERR will be set to one (1) and the program will 
stop upon returning to RDUSER. A MATRD flow chart and the subprogram listing follow. 
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non on o o o o o 


C 

C 

C 

c 

c 

c 

c 


SUBROUTINE MATRD ( NAM 1 « NAM2 * I ERR . I OPT ) 
THIS SUBROUTINE READS MATRIX DATA BLOCKS 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DIMENSION STATEMENTS MUST BE CHANGED 
THE SIZE OF THE MATRICES BEING READ 


IN ACCORDANCE WITH 


DOUBLE PRECISION DVALUOO) 

DIMENSION WRD I N ( 1 00 ) . IRDINC100). VAL(IOO) 

EQUIVALENCE <WRD!N< 1 )* I RDIN ( 1 > ) 

LOGICAL INTGER 

DATA ALP I • SQ 1 . SQ2 ♦ REC 1 .REC2.SYM1 * SYM2 tRE • COM 1 .C0M2/1HI *4H 
1 E * 4HRECT t AHNGLE » 4HSYME . 4HTR I C 1 4HREAL « 4HC0MP . 4HLEX / 

DATA BLANK/4H / 

DATA SING1 * S I NG2 * DOUBl .D0UB2/4H S IN* 4HGLE * 4H DOU t 4HBLE / 


MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

MAT 

SQ U . 4 HA RM A T 
MAT 
MAT 
MAT 


MAT 

******************************************************************* AT 


10 

11 

12 

13 

14 

15 

16 
17 
20 
30 
40 
50 
60 
70 
80 
90 

100 

110 





FORMAT 

FOR A MATRIX DATA BLOCK 



MAT 

1 20 










MAT 

130 

RECORD NO. 

* WORD 

NO 

. * DESCRIPTION 

* 

TYPE 


MAT 

140 

1 

* 

1-2 


* MATRIX LABEL 

* 

ALPHA-NUMERIC 

* 

MAT 

150 

2 

* 


1 


* EOR (-1) 

* 

INTEGER 

♦ 

MAT 

160 

3 

* 


1 


*N0. WORDS NEXT REC. 

* 

INTEGER 

* 

MAT 

170 

4 

* 


1 


* ( TRAILER) GINO NAME 

* 

INTEGER 

* 

MAT 

180 

4 

* 


2 


4 NO. OF COLUMNS 

* 

INTEGER 

♦ 

MAT 

190 

4 

* 


3 


* NO. OF ROWS 

* 

INTEGER 

* 

MAT 

200 

4 

* 


4 


* FORM OF MATRIX 

* 

integer 

♦ 

MAT 

210 

4 

* 


5 


* TYPE OF MATRIX 

* 

INTEGER 

* 

MAT 

220 

4 

* 


6 


*N0. NONZERO TERMS 

♦ 

INTEGER 

♦ 

MAT 

230 

4 

* 


7 


* PER CENT FULLNESS 

♦ 

INTEGER 

* 

MAT 

240 

5 

* 


1 


* EOR (-2) 

♦ 

INTEGER 

♦ 

mat 

250 

6 

* 


1 


*N0. WORDS NEXT REC. 

♦ 

INTEGER 

* 

MAT 

260 

7 

* 


1-2 


* MATRIX LABEL 

♦ 

ALPHA-NUMERIC 

* 

MAT 

270 

8 

* 


1 


* EOR (-3) 

* 

INTEGER 

♦ 

MAT 

280 

9 

* 


I 


♦NO • WORDS NEXT REC. 

* 

integer 

* 

MAT 

290 

10 

* 


1 


♦FIRST NON-ZERO COL. 

* 

INTEGER 

* 

MAT 

300 

10 

* 


2 


* PRECISION 

* 

INTEGER 

* 

MAT 

310 

10 

* 


3 


♦ NOT USED 

* 

INTEGER 

♦ 

MAT 

320 

10 

* 


4 


♦ NOT USED 

* 

INTEGER 

♦ 

MAT 

330 

10 

* 


5 


* NOT USED 

* 

INTEGER 

♦ 

MAT 

340 

10 

* 

6 -NO 

• WORDS 

— 1 * ELEMENTS OF ROW 

♦ 

REAL 

♦ 

mat 

350 

10 

* 

6 -NO 

•WORDS- 

1*P0 INTERS ARE IMBEDDED* 

INTEGER 

* 

MAT 

360 

10 

* 

NO 

• WORDS 

* END OF COLUMN 

* 


♦ 

MAT 

370 

RECORDS 8.9 

'• AND 10 

ARE REPEATED FOR EACH COLUMN 

OF THE MATRIX 


MAT 

380 

WITH 

THE RECORD 

CORRESPONDING to RECORD 8 decreasing by one 


MAT 

390 

EACH 

TIME. 

THIS 

IS 

STOPPED! WHEN A ZERO (0) IS ENCOUNTERED IN 


MAT 

400 

THE 

RECORD 

containing 

THE NUMBER OF WORDS IN 

THE 

NEXT RECORD. 


MAT 

41 0 








i 


MAT 

420 


10 


******* *********************************************#*************MAT 430 

MAT 440 

I CODE= 1 MAT 450 

WRITE (6.380) MAT 460 

I ERR=0 MAT 470 

MAT 480 

READ TAPE MAT 490 

MAT 500 

READ (9) ICNT MAT 510 

IF ( I CNT • EQ. 0 ) GO TO 360 MAT 51 1 

READ <9) ( IRDIN( I ) ♦ 1=1 t I CNT ) MAT 530 
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c 

C CHECK FOR CORRECT MATRIX LABEL 

C 

C WRITE INFORMATION ON BCD TAPE 

IF (NAM1 .EQ. IRDIN <1 ) .AND. NAME. EQ. IRDIN(2) ) GO TO 50 
40 READ (9) JCNT 
READ (9) ICNT 
IF (ICNT.EQ.O) GO TO 10 
READ (9) ( IRDIN ( I) . I =1 , ICNT) 

GO TO 40 

50 IF (ICNT.EQ.2) WRITE (8.390) ( I RD I N (I ) , I * 1 , I CNT ) . I CODE 
READ (9) JCNT 
READ (9) ICNT 
C 

READ INFORMATION FROM TRAILER 

READ (9) (IRDIN! I ), 1*1 * ICNT) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8.400) ( IRDIN ( I ). 1*1 • ICNT > 

I COL* I RD IN ( 2 ) 

I ROW* I RDI N ( 3 ) 

I FORM® I RD I N ( 4 ) 

I TYPE* I RD I N ( 5 ) 

NUMNZ * I RD I N ( 6 ) 

FULL*FLOAT( IRDIN(T) )/I00. 

I CR0W=2* I ROW 

PRINT INFORMATION FROM TRAILER 

WRITE (6.410) NAMI .NAM2. ( IRDIN ( I J) . I J=1 . ICNT ) 

IF (IFORM.EQ.l) FORM 1 = SQ1 
IF (IFORM.EQ.l) F0RM2*SQ2 
IF ( I FORM .EQ .2 ) F0RM1*REC1 
IF (IFORM.EQ.2) F0RM2*REC2 
IF ( I FORM .EQ .6 ) F0RM1*SYM1 
IF ( I FORM .EQ .6 ) F0RM2*SYM2 
IF < I TYPE.EQ . 1 . OR • I TyPE.EQ .2 ) TYPE1=RE 
IF ( I TYPE.EQ. 1 .OR . I TyPE.EQ .2 ) TYPE 2= BLANK 
IF ( I TyPE.EQ *3. OR . I TyPE.EQ .4 ) TYPE1=C0M1 
IF ( I TYPE.EQ .3. OR • I TyPE.EQ .4 ) TYPE2=C0M2 

PRINT MATRIX HEADER 

WRITE (6.420) NAM 1 ,NAM2 . TYPE1 . TYPE2 . ICOL . 1 ROW . FORM 1 .F0RM2 

READ TAPE 

60 READ (9) JCNT 
READ (9) ICNT 

READ < 9 ) ( IRDIN ( I ). 1*1 . ICNT ) 

CHECK FOR CORRECT MATRIX LABEL 
C 

IF (NAMI .EQ. IRDIN(1 > . AND. I RDI N ( 2 ) .EQ.NAM2 ) GO TO 65 
GO TO 360 

65 WRITE (8,460) ICNT 

WRITE (8.550) ( I RDI N ( I ) . I = 1 . I CNT > 

C 

C read tape 

c 


MAT 550 
MAT 560 
MAT 570 
MAT 580 
MAT 590 
MAT 600 
MAT 610 
MAT 620 
MAT 630 
MAT 640 
MAT 650 
MAT 660 
MAT 670 
MAT 680 
MAT 690 
MAT 700 
MAT 710 
MAT 720 
MAT 730 
MAT 740 
MAT 750 
MAT 760 
MAT 770 
MAT 780 
MAT 790 
MAT 800 
MAT 810 
MAT 820 
MAT 830 
MAT 840 
MAT 850 
MAT 860 
MAT 870 
MAT 880 
MAT 890 
MAT 900 
MAT 910 
MAT 920 
MAT 921 
MAT 922 
MAT 923 
MAT 924 
MAT 930 
MAT 940 
MAT 950 
MAT 960 
MAT 970 
MAT 980 
MAT 990 
MATl" 00 
MAT1 01 0 
MAT 1 020 
MAT 1 030 
MAT 1040 
MAT 1 050 
MAT 1 060 
MATl 070 
MAT 1 071 
MAT 1072 
MAT 1 080 
MAT 1 090 
MATl 100 


10 



non 


70 READ (9 > JCNT 
NULC0L=-JCNT-2 
READ <9) ICNT 

TEST FOR ZERO TO END THE READING OF THE MATRIX 

IF (ICNT.EQ.O) GO TO 370 
C 

C TEST FOR NULL COLUMN 

C 

IF ( ICNT.NE. 1 ) GO TO 80 
READ (9) INUMB 
WRITE (8*460) ICNT 
WRITE (8.460) INUMB 
IF (IOPT.EQ.I) PRINT 430* NULCOL 
GO TO 70 
C 

C READ RECORDS OF INFORMATION 

C 

80 READ (9) ( I RDIN ( I ) * I = ! * 5 ) * (WRDIN ( I ) * 1=6 ♦ I CNT ) 

NZROW= IRDIN( 1 ) 

I PREC = I RD I N ( 2 ) 

K = 0 
IK=-3 
C 

C STORE ZEROS IN APPROP1ATE LOCATIONS 

C 

IF (NZR0W.EQ.1) GO TO 100 
C 

C STORE ZEROS IN THE BEGINNING OF THE COLUMN 

C 

NZR0W=NZR0W-1 

K*-l 

DO 90 1=1 *NZROW 
I K= I K+4 
K = K+2 

IF (IPREC.EQ.l) VAL ( I ) =0 « 

IF (IPREC.EQ.2) DVAL (I ) =0*000 
IF (IPREC.EQ.3) VAL(K)=0» 

IF (IPREC.EQ.3) VAL ( K+ 1 ) =0 • 

IF (IPREC.EQ.4) DVAL (K ) =0 .000 
IF (IPREC.EQ.4) DVAL (K+l jsO.ODO 
C 

C STORE DATA FOR BCD TAPE 

C 

IF (IPREC.EQ.2) VAL(K)=0. 

IF (IPREC.EQ.2) VAL ( K+ 1 )=0. 

IF (IPREC.EQ.4) VAL(IK>=0. 

IF (IPREC.EQ.4) VAL(IK+1)=0. 

IF (IPREC.EQ.4) VAL ( IK +2 ) =0. 

IF (IPREC.EQ.4) VAL ( I K+3 ) =0 • 

90 CONT I NUE 

IF (IPREC.EQ.2) K=K+1 
NZROW=NZROW+l 
100 I=NZR0W-1 

IF ( IPREC.EQ.3. OR. IPREC.EQ.4 ) I =2* (NZROW- 1 ) - 1 
IF (IPREC.EQ.2) GO TO 150 
IF (IPREC.EQ.3) GO TO 200 
IF (IPREC.EQ.4) GO TO 250 


MATH 12 
MATl 120 
MAT1 130 
MATl 140 
MAT1 150 
MAT1 160 
MAT1 170 
MAT1 180 
MATl 190 
MAT 1 200 
MAT 1210 
MAT 1 220 
MAT 1 230 
MAT 1 240 
MAT 1 250 
MAT 1260 
MAT 1 270 
MAT 1 280 
MAT 1 290 
MAT 1300 
MAT 1 320 
MAT 1330 
MAT 1340 
MAT 1350 
MAT 1360 
MAT 1370 
MAT 1380 
MAT 1390 
MAT 1 400 
MAT 141 0 
MAT 1420 
MAT 1430 
MAT 1440 
MAT 1 450 
MAT 1 460 
MAT 1 470 
MAT 1480 
MAT 1 490 
MAT 1 500 
MAT 1510 
MAT 1 520 
MATl 530 
MAT 1 540 
MAT 1550 
MAT 1560 
MAT 1570 
MAT 1580 
MAT 1590 
MAT 1 600 
MAT 1610 
MAT 1620 
MAT 1 630 
MAT 1640 
MAT 1650 
MAT 1660 
MAT 1670 
MAT 1680 
MAT 1690 
MAT 1700 
MAT 1710 
MAT 1 720 
MAT 1 730 


C 

C 

C 


STORE SINGLE PRECISION MATRIX ELEMENTS 


a 



o o o n no on o o no o o o o on non 


I WRD*5 

J«NZR0W+ICNT-7 
1 10 1 = 1+1 

I WRD®! WRD+1 

CHECK POINTER FOR IMBEDDED ZEROS 

IF < 1 NIGER ( WRD I N { I WRD ) ) ) GO TO 120 
GO TO 140 

20 JJ=f RDIN< I WRD ) — 1 
1 = 1-1 
30 1=1+1 

STORE IMBEODED ZEROS 
VAL( 1 1=0. 

IF CI.NE.JJ) GO TO 130 
1*1 + 1 

I WRD* I WR0+ 1 

40 VAU ( I ) =WR0IN < IWRD ) 

KCNT* I WR0+NZR0W-6 
IF (KCNT »LT • J ) GO TO 110 
GO TO 300 

STORE DOUBLE PRECISION MATRIX ELEMENTS 

150 IWRD=A 
IK=K-1 

J=NZR0W+ I CNT-6 
160 1=1+1 

I WRD* I WRD+2 

CHECK IMBEDDED ZEROS 

IF ( I NTGER < WRD I N < IWRD) ) > GO TO 170 
GO TO 190 

170 J J= I RO I N ( I WRD ) — 1 
1 = 1-1 

180 1 = 1+1 

STORE IMBEDDED ZEROS 

i ,i 

OVAL ( I )f 0.0D0 

STORE DATA FOR BCD TAPE 

IK= IK+2 
VAl < IK ) *0. 

VAL(IK+1)=0. 

IF ( I .NE. JJ ) GO TO 180 
I » 1 + 1 

IWRD=IWRD+t 

190 DVAl I I ) *DBLE < WRDI N ( I WRD ) ) 

STORE DATA FOR BCD TAPE 
IK= IK+2 

VAL ( I K ) *WRD I N ( I WRD ) 

VAL < IK+I ) *WRD IN ( IWRD+1 1 
KCNT* I WRD+NZROW-4 
IF (KCNT.LT.J1 GO TO 160 
GO TO 300 


MAT 1 740 
MAT 1 750 
MAT 1 760 
MAT 1770 
MAT 1780 
MAT 1790 
MAT 1800 
MAT 1810 
MAT 181 1 
MAT 1 820 
MAT 1 830 
MAT 1 840 
MAT 1 850 
MAT 1860 
MAT 1870 
MAT 1 880 
MAT 1890 
MAT 1900 
MAT 1 91 0 
MAT 1 920 
MAT 1 930 
MAT 1 940 
MAT 1950 
MAT 1 960 
MAT 1970 
MAT 1980 
MAT 1 990 
MAT2000 
MAT201 0 
MAT2020 
MAT2030 
MAT2 _40 
MAT2050 
MAT2060 
MAT2070 
MAT2071 
MAT2080 
MAT2090 
MAT2100 
MAT21 10 
MAT2120 
MAT2I30 
MAT2140 
MAT21 50 
MAT2160 
MAT2170 
MAT2180 
MAT2190 
MAT2200 
MAT2210 
MAT2220 
MAT2230 
MAT2240 
MAT2250 
MAT2260 
MAT2270 
MAT2280 
MAT2290 
MAT2300 
MAT2310 
MAT2320 
MAT2330 
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c 


MAT2340 

c 

STORE COMPLEX MATRIX ELEMENTS 

MAT2350 

c 


MAT2360 

200 

I WRD=4 

MAT2370 


J=NZR0W+ICNT-6 

MAT2380 

2 i 0 

1=1+2 

MAT2390 


I WRD = I WRD+2 

MAT2400 

C 


MAT2410 

c 

CHECK IMBEDDED ZEROS 

MAT2420 

c 


MAT2430 


IF < I NT GER ( WRD I N ( I WRD ) ) ) GO TO 220 

MAT2440 


GO TO 240 

MAT244I 

220 

IF < IRDINI IWRD) .EQ.O > GO TO 240 

MAT2450 


JJ=2*< I RD I N ( I WRD ) -1 )-l 

MAT2460 


1 = 1-2 

MAT2470 

C 


MAT2480 

C 

STORE IMBEDDED ZEROS 

MAT2490 

c 


MAT2500 

230 

1 = 1+2 

MAT2510 


VAL ( I ) *0 • 

MAT2520 


VAL ( 1+1 )=0» 

MAT2530 


IF <I»NE«JJ) GO TO 230 

MAT2540 


1 = 1+2 

MAT2550 


I WRD = I WRD+ 1 

MAT2560 

240 

VAL < I )*WRDIN< IWRD) 

MAT 2570 


VAl ( I + 1 ) =WRD IN( 1 WRD+ I ) 

MAT2580 


KCNT = I WRD+NZROW-4 

MAT2590 


IF IKCNT.LT.J) GO TO 210 

MAT2600 


GO TO 300 

MAT261 0 

C 


MAT2620 

c 

STORE COMPLEX DOUBLE PRECISION MATRIX ELEMENTS MAT2630 

c 


MAT 2640 

250 

IWRD=2 

MAT2650 


J=NZROW+ I CNT— 4 

MAT 2660 

260 

1 = 1+2 

MAT2670 


I WRD= I WRD+4 

MAT 2680 

C 


MAT2690 

C 

CHECK POINTER FOR IMBEDDED ZEROS 

MAT2700 

C 


MAT 271 0 


IF ( INTGER < WRD I N ( IWRD) ) ) GO TO 270 

MAT2720 


GO TO 290 

MAT 2721 

270 

IF ( IRDINI IWRD) .EQ.O ) GO TO 290 

MAT2730 


JJ=2*( IRDINI I WRD ) -1 >-l 

MAT2740 


1 = 1-2 

MAT2750 

C 


MAT 2760 

c 

STORE IMBEDDED ZEROS 

MAT 2770 

c 


MAT2780 

280 

1 = 1+2 

MAT 2790 


DVAL( I )=0.0D0 

MAT.2-.800 


DVAl< 1+1 )=0.000 

MAT2810 

C 


MAT2820 

C 

STORE DATA FOR BCD TAPE 

MAT2830 

C 


MAT2840 


I K= I K+4 

MAT2850 


VAL< IK)=0» 

MAT2860 


VAL < IK+1 ) =0 • 

MAT2870 


VAL< IK+2)*0. 

MAT2880 


VAL< I K+3 ) =0» 

MAT 2890 


IF (I.NE.JJ) GO TO 280 

MAT2900 


1 = 1+2 

MAT2910 


I WRD= I WRD+ 1 

MAT2920 

290 

DVAl( I >=DBLE(WRDIN( IWRD) ) 

MAT2930 
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c 

c 

c 


300 


C 

C 

C 


C 

c 

c 


310 

320 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


DVAL ( 1 + 1 ) = DBLE ( WRD I N ( I WRD+2 ) ) 

STORE DATA FOR BCD TAPE 
I K= IK+4 

VAL( IK)=WRDIN( I WRD) 

VAl(1K+1 )=WRDIN( IWRD+1 ) 

VAL( IK+2)=WRDIN( I WRD+2 ) 

VAL ( IK+3 ) -WRD1N ( IWRD+3) 

KCNT= 1 WRD+NZROW 
IF {KCNT.lt. J) GO TO 260 

IF ( < IPREC.EQ.l .OR. I PREC.EQ.2 > .AND. IROW.EQ. I) GO TO 320 
IF (IPREC.EQ.l. OR. IPREC.EQ.2) K=(I+l)/2 

IF (( IPREC.EQ. 3. OR. IPREC.EQ. 4 ) .AND. IROW.EQ.K ) GO TO 320 
IF ( IPREC.EQ. 3. OR. IPREC.EQ.4 ) K= I 


STORE ZEROS AT THE END OF THE COLUMN 


J= 1 + 1 

DO 310 I B J« I ROW 
IF (IPREC.EQ. 2) 
IF (IPREC.EQ.4) 
K=K+2 

IF (IPREC.EQ.l) 
IF (IPREC.EQ. 2) 
IF (IPREC.EQ. 3) 
IF (IPREC.EQ. 3) 
IF (IPREC.EQ.4) 
IF (IPREC.EQ.4) 


I K= I K+2 
I K= I K+4 

VAL ( I ) =0 « 

DVAL ( I ) =0 . 0D0 
VAL (K > B 0. 

VAL (K+l )®0. 
DVAL <K ) *0 • 0D0 


DVAL (K+l ) =0 . 0D0 
STORE DATA FOR BCD TAPE 


IF (IPREC.EQ. 2) VAL (10*0. 

IF (IPREC.EQ. 2) VAL ( IK+1 ) «0« 

IF (IPREC.EQ.4) VAL(IK)=0. 

IF (IPREC.EQ.4) VAL ( IK + 1 )=0. 

IF (IPREC.EQ.4) VAL ( IK+2 )-0. 

IF (IPREC.EQ.4) VAL ( IK+3 )*0. 

CONTINUE 

IF (IOPT.EQ.l) PRINT 440, NULCOL 

IF (IPREC.EQ. 2) GO TO 330 

IF (IPREC.EQ. 3) GO TO 340 

IF (IPREC.EQ.4) GO TO 350 

PRINT SINGLE PRECISION NUMBERS 

IF (IOPT.EQ.l) WRITE (6,450 ) ( VAL ( I ) , I = 1 , I ROW ) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,460) I ROW 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,470) ( IRDIN ( I ) » 1=1 ,5) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,480) { VAL ( I ) , I B 1 . I ROW ) 

GO TO 70 

PRINT DOUBLE PRECISION NUMBERS 


MAT2940 

MAT 2950 

MAT2960 

MAT2970 

MAT2980 

MAT2990 

MAT3000 

MAT301 0 

MAT3020 

MAT3030 

MAT3040 

MAT3050 

MAT3060 

MAT3070 

MAT3080 

MAT3090 

MAT3100 

MAT31 10 

MAT3120 

MAT3130 

MAT3140 

MAT3150 

MAT3160 

MAT3170 

MAT3180 

MAT3190 

MAT3200 

MAT3210 

MAT3220 

MAT3230 

MAT3240 

MAT3250 

MAT3260 

MAT3270 

MAT3280 

MAT3290 

MAT3300 

MAT331 0 

MAT3320 

MAT3330 

MAT3340 

MAT3350 

MAT3360 

MAT3370 

MAT3380 

MAT3390 

MAT3400 

MAT3410 

MAT3420 

MAT3430 

MAT3440 

MAT3450 

MAT3460 

MAT3470 

MAT3480 

MAT3490 

MAT3500 

MAT3510 

MAT3520 

MAT3530 

M AT3540 

MAT3550 

MAT3560 
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330 


360 


IF (IOPT.EQ.l) WRITE (6,490) ( OVAL ( I ) , I = 1 , I ROW ) 

I ROW 1=1 R0W*2 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,460) IROW1 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,470) ( I RD IN (I), 1=1,5) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,480) ( VA|_ (I), 1 = 1, I ROWl > 

GO TO 70 

PRINT COMPLEX SINGLE PRECISION NUMBERS 

IF (IOPT.EQ.l) WRITE (6,500) ( VAL ( I ) «VAL( 1+1 ) «ALPI , 1=1 , I CROW, 2) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,460) ICROW 

WRITE INFORMATION ON BCD TAPE 


MAT3570 
MAT3580 
MAT3590 
MAT3600 
MAT 361 0 
MAT3620 
MAT3630 
MAT3640 
MAT3650 
MAT3660 
MAT3670 
MAT3680 
MAT3690 
MAT3700 
MAT 37 10 
MAT3720 
MAT3730 
MAT3740 
MAT3750 
MAT3760 
MAT3770 
MAT37B0 
MAT3790 
MAT3800 
MAT 381 0 


WRITE (8,470) ( IRDl N (I ) , I =1 ,5) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,480) ( VAL ( I ) * I = 1 , I CROW ) 

GO TO 70 

PRINT COMPLEX DOUBLE PRECISION NUMBERS 

IF (IOPT.EQ.l) WRITE (6,510) < DVAL ( I ) *DVAl ( I + 1 > » ALP I 

ICROWl=ICROW*2 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,460) ICROW1 

WRITE INFORMATION ON BCD TAPE 

WRITE (8.470) ( IRDIN( I ) , 1=1 ,5) 

WRITE INFORMATION ON BCD TAPE 

WRITE (8,480) (VAL ( I ) , I * 1 , I CROWl ) 

GO TO 70 

RETURN WITH ERROR MESSAGE 

I ERR = I 
RETURN 

PRINT MATRIX INFORMATION 


MAT3820 
MAT3830 
MAT3840 
MAT3850 
MAT3860 
MAT3870 
MAT3880 
MAT3890 
MAT39Q0 
MAT3910 
1=1, ICROW, 2) MAT3920 
MAT3930 
MAT3940 
MAT3950 
M AT3960 
MAT3970 
MAT3980 
MAT3990 
M AT4 ~ 00 
MAT4010 
MAT4020 
MAT4030 
MAT 4040 
MAT4050 
MAT4060 
MAT 4061 
MAT4062 
MAT 4063 
MAT4070 
MAT 4080 
MAT4090 
MAT 41 00 
MAT4110 
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370 WRITE (6.520) NUMNZ 
WRITE (6.530) FULL 

IF ( IPREC.EQ. 1 .OR. I PREC.EQ.3 ) WRITE (6.540 ) SING1.SING2 
IF (IPREC.EQ. 2. OR. IPREC.EQ. 4) WRITE (6.540) 00UB1.D0UB2 

WRITE INFORMATION ON BCD TAPE 
WRITE (8.460) ICNT 
RETURN 
C 

380 FORMAT (1H1) 

390 FORMAT ( 50X. 2 ( A4 . 6X ) . I 1 0 ) 

400 FORMAT (10X.7I10) 

410 FORMAT ( //28H THE TRAILER FOR THE MATRIX .2A4.4H IS .717) 

420 FORMAT (//8H MATRIX .2A4.6H IS A ,2A4. 14. 10H COLUMN X .14. 5H ROW 
12A4.8H MATRIX /) 

430 FORMAT (/8H COLUMN .I6.8H IS NULL/) 

440 FORMAT (/8H COLUMN .16) 

450 FORMAT ( 6E20 . 1 2 . 1 X/ > 

460 FORMAT (70X.I10) 

470 FORMAT (30X.5UO) 

480 FORMAT (4E20.12) 

490 FORMAT (4D30.22. 1 X/ ) 

500 FORMAT ( 3 ( E20 . 1 2 » E20 .1 2 . A 1 ) . 1 X/ ) 

510 FORMAT ( 2 ( D30 .22 . D30 .22 , A 1 ) . 1 X/ ) 

520 FORMAT (/54H THE NUMBER OF NON-ZERO WORDS IN THE LONGEST RECORD 
1.16) 

530 FORMAT (/31H THE DENSITY OF THIS MATRIX IS .F8.2.9H PERCENT /) 
540 FORMAT ( / 1 5H THIS MATRIX I S . 2A4 , 9HPREC I S I ON ) 

550 FORMAT (2 ( A4 »6X ) . (6110)) 

END 


MAT 4120 
MAT 4 1 30 
MAT4140 
MAT 4 150 
MAT4160 
MAT 4170 
MAT4180 
MAT41 <?0 
MAT 4200 
MAT 42 1 0 
MAT 4220 
MAT 4230 
MAT 4240 
. MAT 4250 
MAT4260 
MAT427C 
MAT 4280 
MAT 4290 
MAT 4300 
MAT 4310 
MAT 4320 
MAT 4330 
MAT 4340 
MAT 4350 
MAT 4360 
MAT 4370 
MAT 4380 
MAT 4390 
MAT4391 
M AT44Q0— 
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Subprogram TABRD 

TABRD is a FORTRAN subprogram. Its parameters are the table name (NAM1, 
NAM2), an error flag (IERR), and a print option (IOPT) for the user to print or not to 
print the elements of the table. The primary tasks of TABRD are to read the table data 
from unit 9 and write it on unit 8 in a suitable format. TABRD converts integers into 
real numbers to facilitate the transfer. Integers are marked with a code so that they are 
easily recognized by the WRTUSER program. A TABRD flow chart and the subprogram 
listing follow. 



.WRITE TABLE/ 
DATA ON 
UNITS 
l 6 and 8/ 


:hange integeh 

ALUES TO REAL| 
MARK THEM 
[WITH CODE 
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on non non 


SUBROUTINE TABRD (NAM I ,NAM2 . IERR . I OP T ) TAB 

TAB 

THIS SUBROUTINE READS TABLE DATA BLOCKS TAB 

TAB 

TAB 

DIMENSION STATEMENTS MUST BE CHANGED IN ACCORDANCE WITH TAB 

the size! of The tables being read tab 

- * tab 

DIMENSION I RD IN (100) » WRD IN ( 1 00 ) .VAL < 1 .00') TAB 

EQUIVALENCE! IRDINU ) .WRDIN< 1 ) ) TAB 

LOGICAL INTGER TAB 

TAB 

******************************************************************tab 


THIS SUBROUTINE READS TABLE DATA BLOCKS 


IN ACCORDANCE WITH 


DIMENSION STATEMENTS MUST BE CHANGED 
THE SIZE OF THE TABLES BEING READ 


DIMENSION IRDIN ( 1 00 ) «WRDIN< 1 00 I ♦ VAL < 100 ) 
EQUI VALENCE ( I RD I N < 1 > ♦ WRD I N ( 1 > ) 

LOGICAL INTGER 


FORMAT FOR 


TABLE data BLOCK 


c 

RECORD 

NO* 

* WORD 

z 

o 

• 

* DESCRIPTION 

* 

TYPE 


TAB 

70 

c 

1 

* 

1-2 


* TABLE LABEL 

* 

ALPHA-NUMERIC 

* 

TAB 

80 

c 

2 

* 

1 


* EOR (-1 ) 

* 

INTEGER 

* 

TAB 

90 

c 

3 

* 

1 


♦NO. WORDS NEXT REC. 

* 

I NTEGER 

* 

TAB 

100 

c 

4 

* 

1 


♦ TRAILER GINO NAME 

* 

INTEGER 

* 

TAB 

no 

c 

4 

* 

2~N 


♦ MI SC. INFORMATION 

* ‘ " 

integer 

* 

TAB 

120 

c 

5 

* 

1 


♦ EOR (-2) 

* 

INTEGER 

* 

TAB 

130 

c 

6 

* 

1 


♦NO. WORDS NEXT REC. 

* 

integer 

* 

TAB 

140 

c 

7 

* 

1-2 


♦ TABLE LABEL 

* 

INTEGER 

* 

TAB 

150 

c 

8 

* 

1 


♦ EOR C -3 ) 

* 

INTEGER 

* 

TAB 

160 

c 

9 

* 

1 


♦NO. WORDS NEXT REC. 

* 

INTEGER 

* 

TAB 

170 

c 

10 

* 

all 


♦ELEMENTS OF FIRST 

* 

I NTEGER » REAL* 

TAB 

180 

c 

10 

* 

all 


♦ RECORD OF TABLE 

* 

integer. real* 

TAB 

190 

c 

RECORDS 

8 ♦ 9 ♦ 

AND 10 

ARE 

REPEATED FOR EACH RECORD 

OF THE TABLE, 


TAB 

200 


WITH THE RECORD CORRESPONDING TO RECORD 8 DECREASING BY ONE EACH TAB 
TIME. THIS IS STOPPED WHEN A ZERO CO) IS ENCOUNTERED IN THE RECORDTAB 
CONTAINING THE NUMBER OF WORDS IN THE NEXT RECORD. TAB 

TAB 

* ********************************************************* ******* * T A B 

TAB 

ICODE=2 TAB 

WRITE (6.100) TAB 

IERRsO TAB 

TAB 

READ TAPE TAB 

TAB 

READ (9) ICNT TAB 

IF (ICNT.EQ.O) GO TO 80 TAB 

READ (9) < IRDINU ), 1*1 . ICNT) TAB 

TAB 

CHECK FOR CORRECT TABLE LABEL TAB 

TAB 

IF (NAM1 .EQ. IRDIN ( 1 ) .AND. NAM2.EQ • I RD IN ( 2 ) ) GO TO 50 TAB 

READ (9 ) JCNT TAB 

READ (9) ICNT TAB 

IF (ICNT.EQ.O) GO TO 10 TAB 

READ (9) (IRDIN ( I ) . I 3 1 , ICNT) TAB 

GO TO 40 TAB 

TAB 

WRITE INFORMATION ON BCD TAPE TAB 

TAB 

IF (ICNT.EQ.2) WRITE (8,110) ( I RD I N ( I ) , I = 1 . I CNT ) , I CODE TAB 

READ (9) JCNT TAB 

READ (9) ICNT TAB 

TAB 

read trailer tab 

TAB 
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READ (9) ( IRDIN ( I ) * I =1 * ICNT ) 

C 

C WRITE INFORMATION ON BCD TAPE 

C 

WRITE (8*120) < IRDIN ( ! 1 * 1 = 1 * ICNT ) 

C 

C PRINT TRAILER 

C 

WRITE (6*130) NAM 1 |NAM2 « < IRD!N( I ) « 1 = 1 * I CNT ) 

C 

c read tape 

C 

read (9) JCNT 
READ ( 9 > I CNT 

READ ( 9 ) ( I RD I N < I ) * I = 1 * I CNT > 

C 

C CHECK FOR CORRECT TABLE LABEL 

C 

IF (NAM1 •EOo IROINU ) * AND#NAM2 • EG * I RD I N < 2 ) ) GO TO 60 
GO TO 80 
60 NCNT = 0 

WRITE (6*140) ICNT 

WRITE (8*180) ( IRDlNfl ) *1 = 1 * I CNT ) 

C 

C READ TAPE 

C 

70 READ (9) JCNT 
READ (9) ICNT 

C 

C TEST FOR ZERO (0) TO END THE READING OF THE TABLE 

C 

IF (ICNT.ECUO) GO TO 90 
NCNT =NCNT+1 
C 

C read the values OF THE TABLE 

C 

READ (9) ( WRD IN < I ) * I = 1 * I CNT ) 

C 

C TEST ELEMENTS TO BE REAL OR INTEGER 

C 

KCNT = 0 

DO 75 1=1* I CNT 

IF ( INTGER (WRDIN ( I ) ) •AND.WRDINU ) •NE*0«0 > GO TO 71 
GO TO 72 

71 KCNT=KCNT+t 
VAL ( KCNT ) *-9999*99 

TEMPORARILY CHANGE INTEGER to real * STORE CODE TO SHOW CHANGE 
KCNT-KCNT+1 

VAL (KCNT )*FLOAT( IRDIN ( I ) ) 

GO TO 75 

72 KCNT® KCNT +1 
VAL (KCNT ) =WRD I N ( I ) 

75 CONTINUE 
C 

C WRITE INFORMATION ON BCD TAPE 

C 

WRITE (0*140) KCNT 

WRITE (8* 150 ) (VAL ( 1 ) *1 = 1 *KCNT > 

C 

C PRINT THE VALUES OF THE TABLE 

C 


TAB 550 
TAB 560 
TAB 570 
TAB 580 
TAB 590 
TAB 600 
TAB 610 
TAB 620 
TAB 630 
TAB 640 
TAB 650 
TAB 660 
TAB 670 
TAB 680 
TAB 690 
TAB 700 
TAB 710 
TAB 720 
TAB 730 
TAB 740 
TAB 750 
TAB 751 
TAB 752 
TAB 760 
TAB 770 
TAB 780 
TAB 790 
TAB 800 
TAB 850 
TAB 860 
TAB 870 
TAB 880 
TAB 890 
TAB 900 
TAB 910 
TAB 920 
TAB 921 
TAB 922 
TAB 923 
TAB 924 
TAB 925 
TAB 926 
TAB 927 
TAB 928 
TAB 929 
TAB 930 
TAB 931 
TAB 932 
TAB 933 
TAB 934 
TAB 935 
TAB 936 
TAB 937 
TAB 938 
TAB 939 
TAB 940 
TAB 950 
TAB 960 
TAB 961 
TAB 970 
TAB 980 
TAB 990 
TAB 1000 
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IF (IOPT.EQ.l) WRITE (6.160) NCNT . ( I RD I N ( I ) . I = 1 . I CNT ) TAB1010 

IF ( IOPT.EQ.l) WRITE <6, 170) NCNT , ( WRDIN ( I ) . I - 1 . I CNT ) T AB 1 020 

GO TO 70 TAB 1030 

80 I E RR * 1 TAB! 040 

TAB1 041 

WRITE INFORMATION ON BCD TAPE TAB 1042 


90 

WRITE 

(8.140) I CNT 


r 

RETURN 



V 

100 

FORMAT 

( 1 H 1 ) 


110 

format 

(50X.2(A4,6X) .110) 


120 

FORMAT 

( 1 OX. 71 10) 


130 

FORMAT 

(//27H THE TRAILER FOR THE TABLE .2A4.4H 

140 

FORMAT 

( 70 X * I 1 0 ) 


150 

FORMAT 

(4E20.12 ) 


160 

FORMAT 

( / 34 H THE INTEGER VALUES 

OF RECORD NO* ♦ 

170 

FORMAT 

I/31H THE REAL VALUES OF 

RECORD NO. *I6« 

180 

FORMAT 

END 

(2 ( A4»6X ) * <61 10 ) ) 



IS /( 13110)) 


I6.4H ARE/. (13110) 
4H ARE/. (6E20.12 ) ) 


T AB 1 04 3 
T AB 1 044 
TAB 1050 
T AB 1 060 
TAB1 070 
T AB 1 080 
TAB 1090 
TAB1 100 
TABU 1 0 
T AB 1 120 
) T AB 1 130 
TAB! 131 
TAB 1132 
T AB 1 140- 
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Subprogram INTGER 


INTGER is a logical function subprogram. It is written in FORTRAN or assembly 
language, depending upon the computer on which it is being used. It is a computer - 
dependent routine. INTGER has one parameter (WORD) which contains a number either 
real or integer. INTGER tests each of the exponent bits for a zero and returns a .TRUE, 
if each bit is a zero (indicating this number is an integer); otherwise it returns a .FALSE, 
(indicating this number is real). INTGER is called by both MATRD and TABRD. The 
INTGER subprogram listings for each of the three NASTRAN operative computers follow. 
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IDENT 


INTGER 


I NT 1 





INT 

2 

* 



INT 

3 

* 

logical 

FUNCTION INTGER (WORD) 

INT 

4 

* 



INT 

5 


RETURNS 

TRUE IF UPPER 12 BITS ARE ZERO 

INT 

6 

* 

FOR USE 

ON THE CD C 6000 SERIES COMPUTER 

INT 

7 

* 



INT 

8 


ENTRY 

INTGER 

INT 

9 


VFD 

42/OL I NTGER « 1 8/2 

INT 

10 

INTGER 

BSS 

1 E/E LINE 

INT 

1 1 


SA) 

B 1 

INT 

12 


MX2 

12 

1 N T 

13 


BX3 

X2*X 1 

INT 

14 


MX6 

60 

INT 

15 


ZR 

X3* INTGER 

INT 

16 


SX6 

BO 

INT 

17 


EQ 

INTGER 

INT 

18 


END 


INT 
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LOGICAL FUNCTION INTGER (WORD) 

C 

C THIS SUBROUTINE TESTS THE UPPER EIGHT BITS OF A WORD FOR ZEROS 

C IF ALL EIGHT BITS ARE ZERO INTGER RETURNS TRUE OTHERWISE FALSE 

C FOR USE ON THE IBM 360*370 COMPUTER 

C 

LOG I CAL *1 TEST 

EGU I VALENCE ( TEST « ACHK ) 

ACHK-WORD 
IF (TEST) 10*20*10 
10 I NT GER= ©FALSE « 

RETURN 

20 I NT GER= » TRUE ® 

RETURN 

END 


I NT 10 
I NT 20 
I NT 30 
I NT 40 
IBM 4 1 
I NT 50 
I NT 60 
I NT 70 
I NT 80 
I NT 90 
INT 100 
I NT 110 
INT 120 
INT 130 
INT 140 




LOGICAL FUNCTION INTGER (WORD) 



INT 

10 

C 






INT 

20 

C 


THIS SUBROUTINE TESTS THE 

UPPER 

NINE BITS OF A WORD FOR 

ZEROS 

INT 

30 

C 


INTGER RETURNS TRUE IF THE 

WORD 

IS AN INTEGER OTHERWISE 

IT RETURNSINT 

40 

c 


FALSE 




INT 

50 

c 


FOR USE ON THE UNIVAC 1100 

SERIES COMPUTER 


INT 

51 

c 






INT 

60 



I TEST=FLD (0,9. WORD ) 




INT 

70 



IF ( ITEST >10,20,10 




INT 

80 


10 

INTGER=.FALSE. 




INT 

90 



RETURN 




INT 

100 


20 

I NTGER= • TRUE • 




INT 

1 10 



RETURN 




INT 

120 



END 




INT* 

1 30 
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RDUSER USAGE 


Control-Card Operation for RDUSER 

RDUSER is executed on different computers by a different set of control cards. The 
following three sets are acceptable for the indicated computer: 

(1) CDC 6000 series (Control Data 6000 series computer systems) 

JOB,... 

REQUEST, TAPES), HY. number ,ROL (Binary input tape) 

REWIND (TAPE 9) 

REQUEST, TAPE8,HY,X. number ,RIL, initials, identification (BCD output tape) 

X (external) is optional depending on where the tape is being sent. 

RUN (S) 

LGO. 

DROPFIL(TAPE8) 

DROPFIL(TAPE9) 

7 8q (END-OF-RECORD) 

PROGRAM RDUSER(INPUT, OUTPUT, TAPE5=INPUT, 

T APE 6=OUTPUT , TAPE8, TAPE9) 

* -(source deck} 


• {card input} 

6 7 fi (END-OF-FILE) 
y 9 

(2) IBM 360-370 series 
//JOB,... 

//A EXEC FORTRAN H 
//TIME = number 
//SYSIN DD * 

• {source deck} 

//B EXEC LINKGO 
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//GO.FT05F001 DD * 


//GO.FT08F001 DD UNIT=7TRACK,V0L=SER=number, 

// LABE L=(,NL) ,DISP=NE W ,DSN=name , 

// DCB=(RECFM=F,LRECL=136,BLKSIZE=136,TRTCH=ET (BCD output tape) 
//GO.FT09F001 DD UNIT=2400-3 ,VOL=SER=number, 

// LABEL=(,NL) ,DISP=OLD,DSN=name, 

// DCB=(RECFM=VBS,LRECL=84,BLKSIZE=944) (Binary input tape) 

/* 

(3) UNIVAC 1100 series 


g RUN,//., 
g ASG,T 
§ REWIND 
g ASG,T 
g REWIND 
g FOR, IS 

9,T,SAVE05 

9 

8,T,SAVE05 

8 

RDUSER, RDUSER 

{Binary input tape} 
{BCD output tape} 

• {source 

deck} 


g MAP, I 
g XQT 

relocatable element, absolute element 
absolute element 

• {data deck} 



Error Messages Output by RDUSER 

Messages from RDUSER indicating an error are defined in the following list: 

LABEL (name) WAS NOT FOUND ON TAPE 

User tape label is missing on the input tape (unit 9). 

THE NAME (name) WAS NOT FOUND ON TAPE 

The matrix or table named is missing on the input tape (unit 9). 

Restrictions in RDUSER 

The dimensions for DVAL, VAL, IRDIN, and VAL in the two subroutines MATRD and 
TABRD must be large enough to accommodate a single column of the largest matrix and 
table being transferred. 
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The tape must be created on a 7-track tape drive when the receiving installation has 
a CDC or a UNIVAC computer. 


Card Input for RDUSER 
The input for RDUSER is as follows: 

The first N cards are read with a free field format (20A4). These cards are used to 
describe the data blocks being passed between computers. The Nth card must have a blank 
in column 77 and END in columns 78 to 80 to stop the reading of the comments (see sample 
input in the next section). The next card is read with the format (2A4) and will have the 
FORTRAN User Tape Label (left justified). The remaining cards will be read with the 
format (2A4,2I2). Columns 1 to 8 will have the data block name (left justified). Column 10 
will have a code for determining whether a data block is a matrix (1) or a table (2). Col- 
umn 12 will contain a code determining the option of printing (1) or not printing (0) the 
elements of the data block. 


Sample Input for RDUSER 
A sample input for RDUSER is as follows: 


THE MATRICES ARE TO TEST SINGLE « DOUBLE 
MATRICES . TO TEST THE ORDER 
ORDER OF HEADERS AND FINALLY 
NULL COLUMNS 


xxxxxxxx 
ZC0L5 1 1 

ZC0L6 1 0 

ZC0L8 1 0 

ZCOL1 1 0 

ZCOL2 1 1 

ZC0L3 1 1 

B 1 0 

X I 1 

EQEXIN 2 I 

GPDT 2 0 


PRECISION t COMPLEX 
OF OUTPUT 2 ' THE 


AND REAL 


END 
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DESCRIPTION OF WRTUSER AND ITS SUBPROGRAMS 


Program WRTUSER 

WRTUSER is a FORTRAN main program that has as its primary function the calling 
of various subprograms for data manipulation. The flow chart and program listing that 
follow show how WRTUSER controls program operation. 
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c 

C THIS PROGRAM CONVERTS A BCD TAPE INTO A BINARY UNFORMATTED TAPE 

C FOR USE AS INPUT INTO NASTRAN MODULE INPUTT2 

C 

. DIMENSION NAME (2 ) . I TRAIL <7 ) . IRDIN(IOO) 

C 

C CNVSRC CONVERTS THE SOURCE FROM ONE COMPUTER INTO A SOURCE FORM 

C READABLE BY A DIFFERENT COMPUTER 

C 

CALL CNVSRC 
C 

C READ IN COUNTER FOR NUMBER OF DATA BLOCKS FOR THIS LABEL 

C MACH* 1 IMPLIES MACHINE FOR OUTPUT TAPE IS CDC 

C MACH=2 IMPLIES MACHINE FOR OUTPUT TAPE IS IBM 

C MACH*3 IMPLIES MACHINE FOR OUTPUT TAPE IS UN I VAC 

C 

READ (5*40) NAM1 «NAM2. ITOT.MACH 
C 

C CALL SUBROUTINES FOR COMMENT AND HEADER INFORMATION 

C 

PRINT 50 
CALL WRTCOM 
PRINT 50 
CALL HEADER 
1=0 

20 I CnT=2 

WRITE (8) I CNT 
READ (9.60) NAME » I CODE 

WRITE THE DATA BLOCK NAME 

WRITE (8) NAME 
ICNT=-1 

WRITE (8) I CNT 
I CNT = 7 

WRITE (8) I CNT 
READ (9.70) I TRAIL 

WRITE THE TRAILER 

WRITE (8) l TRAIL 
ICNT=-2 

WRITE (8) I CNT 
READ (9.90) I CNT 
WRITE (8) I CNT 

READ AND WRITE THE DATA BLOCK NAME 

READ (9.100) ( IRDIN(L).L=1 . ICNT ) 

WRITE (8) ( IRDIN(L) »L=1 . ICNT ) 

CALL A SUBROUTINE TO WRITE THE MATRICES 
C 

IF (ICODE.EQ.l) CALL M ATWRT ( MACH ) 

C 

C CALL A SUBROUTINE TO WRITE THE TABLES 

C 

IF (ICODE.EQ.2) CALL TABWRT 
1 = 1 + 1 
C 

C CHECK FOR THE END OF THE DATA BLOCKS FOR THAT HEADING 

C 


WRT 1 1 
WRT 1 2 
WRT 13 
WRT 14 
WRT 20 
WRT 21 
WRT 22 
WRT 23 
WRT 24 
WRT 25 
WRT 40 
WRT 50 
WRT 51 
WRT 52 
WRT 53 
WRT 60 
WRT 70 
WRT 90 
WRT 100 
WRT 110 
WRT 120 
WRT 130 
WRT 140 
WRT 150 
WRT 160 
WRT 170 
WRT 1 80 
WRT 190 
WRT 200 
WRT 210 
WRT 220 
WRT 230 
WRT 240 
WRT 250 
WRT 260 
WRT 270 
WRT 280 
WRT 290 
WRT 300 
WRT 310 
WRT 320 
WRT 330 
WRT 340 
WRT 350 
WRT 360 
WRT 370 
WRT 380 
WRT 390 
WRT 400 
WRT 401 
WRT 410 
WRT 420 
WRT 430 
WRT 440 
WRT 450 
WRT 460 
WRT 470 
WRT 480 
WRT 490 
WRT 500 
WRT 510 
WRT 520 
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40 

50 

60 

70 

80 

90 

100 


PRINT 80 « NAME 
IF (l.LT.ITOT) GO TO 20 
END FILE 8 
STOP 

FORMAT (2A4 *215) 

FORMAT ( 1 H 1 ) 

FORMAT ( 50X « 2 ( A4 « 6X ) « I 1 0 ) 

FORMAT ( 1 OX » 7 1 1 0 ) 

FORMAT (1H0.12H DATA BLOCK .2A4.25H HAS BEEN WRITTEN ON TAPE) 
FORMAT ( 70X ♦ 110) 

FORMAT <2( A4«6X ) « <61 10 ) ) 

END 


WRT 530 
WRT 540 
WRT 560 
WRT 570 
WRT 580 
WRT 590 
WRT 600 
WRT 610 
WRT 620 
WRT 630 
TAB 631 
TAB 632 
WRT 640- 
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Subprogram CNVSRC 

CNVSRC is a FORTRAN subprogram. The primary job of CNVSRC is to convert 
the source of one computer into a source form readable by a computer of a different type. 
The CNVSRC subprogram listing follows. 


SUBROUTINE CNVSRC 


CONV 10 


10 


20 


30 

40 

50 

60 


THIS PROGRAM CONVERTS THE SOURCE FROM 
SOURCE FOR ANOTHER COMPUTER (MACH2) 

DIMENSION I DATA (20) 

DATA ISTOP. IT/4HST0P.4H IT / 

K = 0 

READ IN COMPUTER TYPES 

READ (5*30) MACH1 ,MACH2 
PRINT 40.MACH1 .MACH2 

READ IN DATA 

READ (3*50) IDATA 
IF { IDATA < 19) .EQ. ISTOP* AND. IDATA <20 ) *EQ 
K*K+ 1 

WRITE ( 9 ♦ 50 ) IDATA 
GO TO 10 
PRINT 60. < 

REWIND 9 
RETURN 

FORMAT (A4.2X.A4) 

FORMAT < 1H1.34H START CONVERTING SOURCE 

1 TO THE .A4.9H COMPUTER /) 

FORMAT < 2 0 A 4 ) 

FORMAT <//24H STOP SOURCE CONVERSION ./ 
1 CONVERTED ) 

END 


CONV 20 

ONE COMPUTER (MACH! ) TO THE CONV 30 

CONV 40 
CONV 50 
CONV 60 
CONV 70 
CONV 80 
CONV 90 
C0NV100 
CONV1 10 
C0NV120 
C0NV130 
CONV 140 
C0NV150 
C0NV160 
C0NV170 

IT) GO TO 20 CON VI 80 

CONV190 
C0NV200 
C0NV210 
C0NV220 
C0NV230 
CONV240 
CONV250 
CONV260 

FROM THE .A4.17H COMPUTER CONV270 

CONV280 
C0NV290 

X.I5.29H RECORDS HAVE 8EENCONV300 

CONV310 
CONV320 


29 



o u u 


Subprogram WRTCOM 


WRTCOM is a FORTRAN subprogram that reads the comments generated by 
RDUSER on unit 9 and prints them out. The WRTCOM subprogram listing follows* 


SUBROUTINE WRTCOM 

THIS SUBROUTINE READS AND PRINTS COMMENTS 

DIMENSION I COM ( 20 ) 

DATA IQUIT/4H END/ 

10 READ ( 9 « 20 ) { I COM (0*1*1 *20 ) 

PRINT 20« ( I COM ( I > • 1*1 ♦ 20 ) 

IF i I COM ( 20 ) eNE * I QU IT) GO TO 10 
RETURN 

20 FORMAT <20A4 > 

END 


COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 


10 
1 1 
12 
13 
20 
30 
40 
50 
60 
70 
80 
90 

i no- 
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Subprogram HEADER 


HEADER is a FORTRAN subprogram. Its purpose is to write the header informa- 
tion on unit 8 in binary form. The header information includes the data, header, and 
user tape label. The HEADER subprogram listing follows. 


SUBROUTINE HEADER 

SUBROUTINE FOR WRITING HEADER INFORMATION 

DIMENSION I DATE <3 ) * HEAD (7 ) » NAME ( 2 ) 
ICNT=3 

WRITE (8) 1 CNT 

READ (9*10) I DATE 

WRITE THE DATE 

WRITE (8) I DATE 
ICNT=7 

WRITE (8) I CNT 
READ (9*20) HEAD 
C 

C WRITE THE HEADER 

C 

WRITE < 8 ) HEAD 
ICNT=2 

WRITE (8) I CNT 
READ (9*30) NAME 
C 

C WRITE THE LABEL 

C 

WRITE (8) NAME 
ICNT=-I 

WRITE (8) I CNT 
ICNT=0 

WRITE (8) I CNT 
RETURN 
C 

10 FORMAT <50X*3I10) 

20 FORMAT ( 1 OX* 7 ( A4 , 6X > ) 

30 FORMAT ( 60X * 2 < A4 * 6X > > 

END 


HED 1 0 
HED 20 
HED 30 
HED 40 
HED 50 
HED 60 
HED 70 
HED 80 
HED 90 
HED 100 
HED 1 I 0 
HED 120 
HED 130 
HED 140 
HED 150 
HED 160 
HED 170 
HED 180 
HED 190 
HED 200 
HED 210 
HED 220 
HED 230 
HED 240 
HED 250 
HED 260 
HED 270 
HED 280 
HED 290 
HED 300 
HED 31 0 
HED 320 
HED 330 
HED 340 
HED 350 
HED 360- 
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Subprogram MATWRT 


MATWRT is a FORTRAN subprogram. Its parameter is the computer (MACH) on 
which the tape output on unit 8 from WRTUSER will be used. The primary task for 
MATWRT is to pack a matrix a column at a time and write it in binary form on unit 8. 

A MATWRT flow chart and a subprogram listing follow. 
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SUBROUTINE MATwRT ( MACH > 

C 

C SUBROUTINE FOR WRITING MATRICES 

C 

c 

C DIMENSION STATEMENT MUST BE ADJUSTED ACCORDING TO SIZE OF MATRIX 

C 

DIMENSION IRDIN (5 ) » WRDIN(IOO), IVAL(200>. VAL(200) 

EQUIVALENCE ( IvAL (I ) »V4L( 1 ) > 

IF (MACH.EQ.l) I STORE® 131071 
IF (MACH.EQ.2) I STORE® 16777215 
IF (MACH.EQ.3) IST0RE=29843125 
JCNT=-2 

10 JCNT= JCNT— 1 

WRITE (8) JCNT 
READ (9.1701 I CNT 

CHECK FOR END OF DATA BLOCK 

IF (ICNT.EQ.O) GO TO 160 
IF (ICNT.NE.l) GO TO 20 
READ (9.170) INUMB 
WRITE (8) I CNT 
WRITE (8) I STORE 
GO TO 10 

20 READ (9.180) ( IRDIN( I 1 , 1=1 ,5 ) 

NZCOL= IRDIN ( 1 > 

I PREC= 1 RD IN ( 2 ) 

READ (9.190) (WRDIN( I ) . 1=1 . ICNT ) 

KCNT=0 
ISW = 0 

I PREC* 1 IMPLIES SINGLE PRECISION REAL MATRIX 
IPREC=2 IMPLIES DOUBLE PRECISION REAL MATRIX 
I PREC=3 IMPLIES SINGLE PRECISION COMPLEX MATRIX 
I PREC=4 IMPLIES DOUBLE PRECISION COMPLEX 

GO TO (30.70.70.110), IPREC 

WRITE SINGLE PRECISION REAL MATRICES 

30 DO 60 I =NZCOL« ICNT 

CHECK FOR IMBEDDED ZEROS 

IF (WRDIN( I ) .EO.O. ) GO TO 50 
IF ( I SW.EQ.O ) GO TO 40 
KCNT=KCNT+1 

STORE POINTER TO NEXT NON-ZERO ELEMENT 

I VAL (KCNT )= I 
40 <CNT=KCNT+1 
ISW=0 

STORE NON-ZERO ELEMENTS 

VAL(KCNT)=WRDIN(I ) 

GO TO 60 


mat 

10 

MAT 

20 

MAT 

30 

MAT 

40 

mat 

41 

MAT 

42 

MAT 

43 

MAT 

50 

MAT 

60 

MAT 

61 

MAT 

62 

MAT 

63 

mat 

70 

MAT 

00 

MAT 

90 

MAT 

100 

MAT 

110 

MAT 

120 

MAT 

130 

MAT 

140 

MAT 

150 

MAT 

160 

MAT 

170 

MAT 

180 

MAT 

190 

MAT 

200 

MAT 

220 

MAT 

230 

MAT 

240 

MAT 

250 

MAT 

260 

MAT 

270 

MAT 

280 

MAT 

290 

MAT 

300 

MAT 

310 

MAT 

320 

MAT 

330 

MAT 

340 

MAT 

350 

MAT 

360 

MAT 

370 

MAT 

380 

MAT 

390 

MAT 

400 

MAT 

410 

MAT 

420 

mat 

430 

mat 

440 

mat 

450 

MAT 

460 

MAT 

470 

MAT 

480 

MAT 

490 

MAT 

500 

MAT 

510 

MAT 

520 

MAT 

530 

MAT 

540 
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50 I SW*1 
60 CONTINUE 

LCNT=KCNT*6 
WRITE (8) LCNT 

WRITE COLUMN OF MATRIX ON TAPE 

WRITE < S) ( IRDINC I ) * I B 1 «5 ) ♦ ( V AL < I ) * I = 1 * KCNT ) ♦ I STORE 
GO TO 10 

WRITE DOUBLE PRECISION REAL MATRICES 
WRITE SINGLE PRECISION COMPLEX MATRICES 

70 I ZCOL = < 2*NZC0L )-1 

DO 100 I = I ZCOL * I CNT *3 

CHECK FOR IMBEDDED ZEROS 

IF < WRDIN ( I ) ®EQ#0 ® * AND#WRD I N ( ! +1 ) ©EQ«0e ) GO TO 90 

IF < I SW®EGU 0 ) GO TO 80 

KCNT=KCNT+1 

STORE POINTER TO NEXT NON-ZERO ELEMENT 


MAT 550 
MAT 560 
MAT 570 
MAT 580 
MAT 590 
MAT 600 
MAT 610 
MAT 630 
MAT 630 
MAT 640 
MAT 650 
MAT 660 
MAT 670 
MAT 680 
MAT 690 
MAT 700 
MAT 710 
MAT 720 
MAT 730 
MAT 740 
MAT 750 
MAT 760 
MAT 770 


80 


90 

100 


110 


120 


IVALCKCNT >=( 1 + 1 )/2 
KCNT =KCNT + 1 

STORE NON-ZERO ELEMENTS 

VAL ( KCNT ) * WRD I N ( I ) 

KCNT=KCNT+1 

VAL ( KCNT )*WRDIN( 1 + 1 ) 

ISW = 0 

GO TO 100 

ISW=i 

CONTINUE 

LCNT=KCNT+6 

WRITE (8) LCNT 

WRITE COLUMN OF MATRIX ON TAPE 


MAT 780 
MAT 790 
MAT 800 
MAT 810 
MAT 830 
MAT 830 
MAT 840 
MAT 850 
MAT 860 
MAT 870 
MAT 880 
MAT 890 
MAT 900 
MAT 910 
MAT 930 
MAT 930 
MAT 940 


WRITE (8 ) ( IRDINC I ) * 1 = 1*5) * CVALI I ) * 1=1 *KCNT ) * I STORE 

GO TO 10 

WRITE DOUBLE PRECISION COMPLEX MATRICES 

IZCOL= C4*NZC0L>-3 
DO 150 I = I ZCOL * 1 CNT « 4 

CHECK FOR IMBEDDED ZEROS 

IF (WRDIN ( I ) *EQ®0®# AND. WRDIN ( 1+1 ) ®EQ»0 ® ©AND© WRDIN ( I+3)©EQ©0 
1 RD INC 1+3 ) «EQ ® 0 « ) GO TO 140 
IF ( I SW ® EQ ® O ) GO TO 120 
K CNT = KCNT +1 

STORE POINTER TO NEXT NON-ZERO ELEMENT 

I VAL (KCNT )* ■( 1+3 )/A 
t SW=0 

DO 130 J=1 *4 
KCNT=KCNT+1 


MAT 950 
MAT 960 
MAT 970 
MAT 980 
MAT 990 
MAT! 00 
MAT! 010 
MAT 1 020 
MAT1 030 
MAT 1040 
MAT! 050 
> • AND* WMAT 1 060 
MAT 1 070 
MAT! 080 
MAT 1 090 
MAT1 100 
MAT1 1 10 
MAT 1 130 
MAT1 130 
MAT1 140 
MAT1 150 
MAT 1 160 
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c 

C STORE NON-ZERO ELEMENTS 

C 

VAL<KCNT)*WRDINC I+J-l ) 

130 CONTINUE 
GO TO 150 
140 I SW= 1 
150 CONTINUE 

LCNT »KCNT +6 
WRITE (8) LCNT 
C 

WRITE COLUMN OF MATRIX ON TAPE 

WRITE ( 8 > ( I RO IN<I)*I*1»5)«< VAL <I)*I*1« KCNT ) * I ST ORE 

GO TO 10 

160 WRITE (8) ICNT 
RETURN 
C 

170 FORMAT (70X.I10) 

180 FORMAT <30X«5Il0) 

190 FORMAT <4E20»12) 

END 


MAT 11 70 
MAT 1 180 
MAT1 190 
MAT 1200 
MAT 1210 
MAT 1220 
MAT 1230 
MAT 1 240 
MAT 1250 
MAT 1260 
MAT 1270 
MAT 1 280 
MAT 1290 
MAT 1300 
MAT 1310 
MAT 1320 
MAT 1 330 
MAT 1340 
MAT 1350 
MAT 1360 
MAT l 370 
MAT 1380- 
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Subprogram TABWRT 


TABWRT is a FORTRAN subprogram that reads a table from unit 9, converts real 
numbers to integers if needed, and writes the table on unit 8 in binary form. A TABWRT 
flow chart and a subprogram listing follow. 



36 




non non 


SUBROUTINE TABWRT 


C 

C SUBROUTINE FOR WRITING TABLES 

C 

c 

C DIMENSION STATEMENT MUST BE ADJUSTED ACCORDING TO SIZE OF TABLE 

C 

DIMENSION IVAL( 100) , VAL (100) ,WRD!N( 100) 

EQUIVALENCE ( I VAL ( 1 > * VAL ( 1 > > 

JCNT *— 3 

10 WRITE (8) JCNT 
READ (9,30) ICNT 

C 

C CHECK FOR END of DATA BLOCK 

C 

IF (ICNT.EQ.O) GO TO 20 
READ (9,40) <WRDIN< I )* 1*1 , ICNT) 

CHANGE REALS BACK INTO INTEGERS IF NEEDED 

KCNT=0 
1=0 

11 1 * 1+1 

IF ( I .GT. ICNT ) GO TO 13 
IF ( WRDI N ( I ) *EQ» -9999 ,99 ) GO TO 12 
KCNT*KCNT+1 
VAL(KCNT)=WRDIN( I ) 

GO TO 11 
12 1*1+1 

KCNT =KCNT + 1 

lVALCKCNT )=IFIX(WRDIN( 1 ) ) 

GO TO 11 

WRITE THE TABLE 

13 WRITE (8) KCNT 

WRITE (B) ( VAL ( I ) * I = 1 , KCNT ) 

JCNT = JCNT— 1 
GO TO 10 

20 WRITE (8) ICNT 
RETURN 
C 

30 FORMAT (70X»I10> 

40 FORMAT (4E20» 12 ) 

END 


TAB 10 
TAB 20 
TAB 30 
TAB 40 
TAB 41 
TAB 42 
TAB 43 
TAB 50 
TAB 51 
TAB 60 
TAB 70 
TAB 80 
TAB 100 
TAB 110 
TAB 120 
TAB 130 
TAB 140 
TAB 141 
TAB 142 
TAB 143 
TAB 144 
TAB 145 
TAB 146 
TAB 147 
TAB 148 
TAB 149 
TAB 150 
TAB 151 
TAB 152 
TAB 153 
TAB 154 
TAB 155 
TAB 150 
TAB 160 
TAB 170 
TAB 179 
TAB 180 
TAB 190 
TAB 200 
TAB 210 
TAB 211 
TAB 220 
TAB 230 
TAB 240 
TAB 250- 
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Control-Card Operation for WRTUSER 

WRTUSER is executed on different computers by a different set of control cards. 
The following three sets are acceptable for the indicated computer: 

(1) CDC 6000 series 
JOB,... 

REQUEST TAPE3,HY,X. number ,ROL (BCD input tape) 

X (external) is optional depending on where the tape was generated 
REWIND ( TAPE 3 ) 

REQUEST(TAPE8),HY. number ,RIL, initials, identification (Binary output tape) 
REWIND ( TAPE 8) 

RUN(S) 

LGO. 

DROPFIL(TAPE3) 

DROPFIL(TAPE8) 

EXIT. 

DROPFIL(TAPE3) 

DROPFIL(TAPE8) 

7 8 g (END-OF-RECORD) 

PROGRAM WRTUSER(INPUT , OUTPUT , TAPE 5=INPUT , TAPE 8 , TAPE 9 , TAP E 3 ) 

• {source deck} 


7 8 g (END-OF-RECORD) 
• {card input} 


6 7 fl (END-OF-FILE) 
y 9 

(2) IBM 360-370 series 
//JOB,... 

//A EXEC FORTRANH 
/ /TIME=number 
//SYSIN DD * 
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• {source deck} 

//B EXEC LINKGO 
//GO.FT05F001 DD * 

• {card input} 

//GO. FT03F001 DD UNIT= 7 TRAC K, VO L=SER=number , L AB E L= ( ,NL) , 

// DISP=OLD,DSN=name,DCB=(RECFM=F,LRECL=136,BLKSIZE=136,TRTCH=ET) 
(BCD input tape) 

//GO.FT08F001 DD UNIT=2400-3,VOL=SER=number,LABEL=(,NL), 

// DISP=NEW ,DSN=name,DCB=(RECFM=VBS,LRECL=84,BLKSIZE=844) 

(Binary output tape) 

//GO.FT09F001 DD UNIT=2 314 ,DISP=NE W , SPACE =(TRK, (10,1)), 

// DCB=(RECFM=FB,LRECL=80,BLKSIZE=960) 

/* 

(3) UNIVAC 1100 series 


8 RUN, // 
% ASG,T 

3,T,SAVE05 

(BCD input tape) 

g REWIND 

3 


e ASG,T 

8,T,SAVE05 

(Binary output tape) 

@ REWIND 

8 


8 ASG,T 

9,T 

(scratch space) 

8 FOR, IS 

WRTUSER, WRTUSER 



• {source deck} 

■8 MAP, I relocatable element, absolute element 

% XQT absolute element 

• {data deck} 


Error Messages Output by WRTUSER 
There are no error messages in the WRTUSER program. 
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Restrictions in WRTUSER 


The dimensions for IVAL, VAL, and WRDIN in the two subroutines TABWRT and 
MATWRT must be large enough to accommodate a single column of the largest matrix 
and table being transferred. 


Card Input for WRTUSER 
The input for WRTUSER is as follows: 

The first card, a comment card, is read with the format (A4,2X,A4) where columns 1 
to 4 have the type of computer (left justified) on which the tape was created and columns 7 
to 10 contain the type of computer (left justified) on which the tape is being converted. The 
second card is read with the format (2A4,2I5), Columns 1 to 8 have the FORTRAN User 
Tape Label (left justified). Columns 9 to 13 have the number of data blocks to be con- 
verted. Column 18 has a code to show that the tape is being converted for use on a 
CDC (1), IBM (2), or UNIVAC (3) computer. 


Sample Input for WRTUSER 
The following is a sample input for WRTUSER: 

IBM CDC 

XXXXXXXX 10 1 


VERIFICATION OF PROGRAMS 

RDUSER and WRTUSER were executed for four* of the nine possibilities shown in 
figure 1 and found to possess the desired qualities lacking in DMI punched cards. Card 
handling for the input to NASTRAN was cut to the minimum. Square, rectangular, and 
symmetric matrices containing single -precision real, single-precision complex, double- 
precision real, and double -precision complex elements were used in the test runs. In 
each case the answers listed on one computer agreed with the answers listed on dissimilar 
and similar computers, which indicated that no precision was lost in the transfer. 

^UNIVAC paths were not tested due to errors in the INPUTT2 and OUTPUT2 
NASTRAN modules on the UNIVAC computer. 
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Figure 1.- Paths of data between computers. 


CONCLUDING REMARKS 

The transfer of data between dissimilar computers by using the two new NASTRAN 
utility programs RDUSER and WRTUSER has been done successfully for the IBM and CDC 
computers. The transfer of data to and from the UNI VAC computer has not been com- 
pletely tested due to errors in the NASTRAN modules OUTPUT2 and INPUTT2. Square, 
rectangular, and symmetric matrices containing single-precision real, single-precision 
complex, double-precision real, and double-precision complex elements were used in the 
test runs. In each case the answers listed on one computer agreed with the answers listed 
on a similar or dissimilar computer, which indicated that no precision was lost in the 
transfer. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 24, 1973. 
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TABLE I.- FORMAT OF HEADER FOR THE NASTRAN USER TAPE 


Record 

Word 

7af 

Description 

1 

1 

I 

Number of words in next record (3) 

2 

1 to 3 

I 

Date 

3 

1 

I 

Number of words in next record (7) 

4 

1 to 7 

A 

Header 

5 

1 

I 

Number of words in next record (2) 

6 

1 to 2 

A 

User tape label 

7 

1 

I 

End of record (-1) 

8 

1 

I 

End of file (0) 

9 

1 

I 

Number of words in next record (2) 


a A = Alphanumeric 
I = Integer 
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TABLE H.- FORMAT OF MATRIX FOR THE NASTRAN USER TAPE 


Record 

Word 

Type 

(a) 

Description 

i 

1 to 2 

A 

Matrix label 

2 

1 

I 

End of record (-1) 

3 

1 

I 

Number of words in next record (7) 

4 



Trailer 


1 

I 

Gina name 


2 

I 

Number of columns 


3 

I 

Number of rows 

' 


4 

I 

Form of matrix 


5 

I 

Type of matrix 


6 

I 

Number of nonzero terms in the longest record 


7 

I 

Percent fullness of matrix 

5 

1 

I 

End of record (-2) 

6 

1 

I 

Number of words in next record (2) 

7 

1 to 2 

A 

Matrix label 

8 

1 

I 

End of record (-3) 

9 

1 

I 

Number of words in next record 

10 



Column of matrix 


1 

I j 

First nonzero row 


2 

I 

Precision of matrix 


3 to 5 

I 

Not used 


6 through 

R,I 

Elements of column. Integer pointers to nonzero elements 


no. words -1 


are imbedded in this record 


No. words 

I 

End of column 

11 

1 

I 

End of record 

12 to N 



Records 9, 10, and 11 are repeated until a zero (0) is found for 




the record containing the number of words in the next record 

N + 1 

1 

I 

Either a zero (0) for an end of file or a two (2) showing the 




number of words in the next record 


a A = Alphanumeric 
I = Integer 
R = Real 


43 




TABLE III.- FORMAT OF TABLE FOR THE NASTRAN USER TAPE 


Record 

Word 

Type 

(a) 

Description 

1 

1 to 2 

A 

Table label 

2 

1 

I 

End of record (-1) 

3 

1 

I 

Number of words in next record (7) 

4 



Trailer 


1 

I 

Gino name 


2 to N 

I 

Miscellaneous information 

5 

1 

I 

End of record (-2) 

6 

1 

I 

Number of words in next record (2) 

7 

1 to 2 

A 

Table label 

8 

1 

I 

End of record (-3) 

9 

1 

I 

Number of words in next record 

10 

All 

I 

Element of first record of table 

11 

1 

I 

End of record 

12 to N 



Records 9, 10, and 11 are repeated until a zero (0) is found for the 




number of words in the next record 

N + 1 

1 

I 

Either a zero (0) for an end of file or a two (2) showing the number 




of words in the next record 


a A = Alphanumeric 
I = Integer 
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